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Abstract 

^r) ' Given a family of rational curves depending on a real parameter, denned by its 

parametric equations, we provide an algorithm to compute a finite partition of the 
parameter space (K, in general) so that the shape of the family stays invariant along 
each element of the partition. So, from this partition the topology types in the family 
C/3 ■ can be determined. The algorithm is based on a geometric interpretation of previous 

c*2 . work ([1]) for the implicit case. However, in our case the algorithm works directly 

with the parametrization of the family, and the implicit equation does not need to 
be computed. Timings comparing the algorithm in the implicit and the parametric 
cases are given; these timings show that the parametric algorithm developed here 
' provides in general better results than the known algorithm for the implicit case. 
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Given a family T of algebraic curves depending on a real parameter, it is 
clear that the shape of the curves in the family may change as the value of 
the parameter changes. In the C.A.G.D. context, one has a good example of 
this phenomenon in the family of offset curves to a given curve (see [3], 
|14j . [T5]). where the parameter is the offsetting distance. Take for example 
the well-known case of the offset to the parabola (see [6]). In this case three 
different shapes, which can be seen in Figure 1 (together with the original 
parabola, in thinner line), arise; here, one may observe that for distances 
d < 1/2, the offsets show no cusps, but they exhibit an isolated point (i.e. 
a geometric extraneous component), and that for d > 1/2, two cusps and a 
self-intersection arise. 
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Fig. 1. Topology types of the offsets to y 2 — x = 0: d < 1/2 (left); d = 1/2 (center); 
d > 1/2 (right) 

. In the middle of each figure, the parabola has been also plotted. 

The problem of computing the different shapes arising in a family of algebraic 
curves depending on a parameter has been considered in [JJ, [12], under dif- 
ferent perspectives but reaching similar results, for the case when the family 
is implicitly given; so, in these papers the family to be analyzed is defined 
as the zero set of a polynomial F(x,y, A), where A is regarded as a real pa- 
rameter. In particular, using the result in [I], the topology types in Figure 1 
can be determined. For this purpose, the strategy is the following. First, one 
computes a finite set A so that in between two consecutive elements of A, the 
topology type of the family does not change; we refer to such a set as a critical 
set; in the case of the offset to the parabola, A = {1/2}. Then, this critical 
set induces a finite partition of the parameter space; in the above example, 
the partition is (0,1/2), {1/2}, (1/2, oo) (notice that the offsetting distance 
is a positive real number). Finally, taking a representative for each element of 
the partition, and applying if necessary standard methods for describing the 
topology of a plane algebraic curve (see [5], [8], [9]), the different shapes in 
the family can be computed. 

In this paper, starting from the results in [T] we address the same problem but 
for the case, specially interesting in the C.A.G.D. context, of a family T of 
rational curves depending on a parameter, defined in parametric form. Thus, 
our input is 

x = u(t, A) 

< 

y = v(t,\) 

v 

where u(t, A), v(t, A) are rational functions in terms of t, A, and A G R is a 
parameter. In these conditions, we consider the computation of a critical set 
of J 7 , but without computing or making use of the implicit equation of the 
family. So, the main result of the paper is an algorithm for carrying out this 
computation. 

The algorithm is advantageous when a rational parametrization of the family 
to be analyzed is available. In fact, in our experimentation we have found 
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several examples (see the comparison table in Subsection I4.3[) of families that 
are hard to study with the result in [1] (i.e. the computations are too heavy in 
implicit form), but which can be analyzed with the algorithm provided here. 

The algorithm is based on a geometrical interpretation of the results in [TJ. 
This interpretation suggests a certain geometrical process to be performed 
in order to compute a critical set. Hence, the question is to carry out that 
process in parametric form. The geometric analysis of the main results in 
[TJ can be found in Section 2. In Section 3, we show how to check different 
hypotheses that we impose on the parametrization, and how to "prepare" the 
family before applying the algorithm; the ideas in this section are illustrated 
in Example [TJ where the offset family to a cardioid is considered. In Section 
4, we provide the full algorithm, we complete the analysis of Example [TJ, and 
we give a comparison table showing the timings of our algorithm, compared 
to those corresponding to the algorithm in [TJ; these timings show that our 
algorithm yields in general better results than the known algorithm for the 
implicit case. Finally, in Section [5] we present the conclusions of our study. 
The parametrizations used for comparing timings are given in Appendix I. 



2 Preliminaries. 

2.1 Known Results for the Implicit Case. 

Let F G R[a;,2/,A]. For all Ao € R such that F(x,y,Xo) is not identically 
we have that F(x,y,X ) = F\ a (x,y) defines an algebraic curve. So, we can 
say that F defines a family T of algebraic curves algebraically depending on 
the parameter A. By Hardt's Semialgebraic Triviality Theorem (see Theorem 
5.46 in [1]), it holds that the number of topology types of T (i.e. the different 
shapes arising in the family as Ao moves in R) must be finite. Therefore, it 
makes sense to consider the problem of computing the topology types arising 
in the family. This problem has been addressed in [TJ, for the case when the 
family is implicitly defined. So, in this subsection we recall the main aspects 
(hypotheses, notation, and results) of this paper. 

More precisely, given F(x, y, A) one can associate an algebraic surface S with 
the family T defined by F by substituting A := z in F; thus, the members 
of T are the level curves of S, i.e. the sections of S with planes normal to 
the z-axis. So, the problem of computing the topology types in T can be re- 
interpreted as the computation of the topology types arising in the family of 
level curves of S. This is exactly the question addressed in [TJ. In order to 
solve this problem, in [TJ the following definition is introduced. 
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Definition 1 Let S be an algebraic surface. We say that A C M is a critical 

set of S, if it is finite and it contains all the z-values where the topology type 
of the level curves of S changes; i.e., if the topology type of the level curves of 
S stays invariant along any interval delineated by two consecutive elements of 



We speak about the critical set of JF, to mean the critical set of the surface 
associated with the family. Notice that if a critical set A = {ai, . . . ,a r } is 
computed, then the parameter space (R, in this case) can be decomposed as 

(—00, ai) U {ai} U (ai, 02) U ■ ■ • U {a r } U (a r , +00) 

Then, taking a representative for each element of the above partition, and 
applying standard methods ([S], [5], for describing the topology of an al- 
gebraic curve, the topology types in the family can be computed. Hence, the 
crucial question is the computation of a critical set. This is the problem ad- 
dressed in p] . In the rest of the subsection, we recall the hypotheses, notation 
(that we will also follow here) and main result of pQ. 

Hypotheses: The hypotheses imposed in [1] on the surface S to be analyzed, 
are: (a) F is square-free, and depends on the variable y; (b) F does not contain 
any factor only depending on the variable z (i.e. S has no component normal 
to the z-axis); (c) the leading coefficient of F w.r.t. the variable y, does not 
depend on x. 

Notation: Given a polynomial G, y/G denotes the square-free part of G, i.e. 
the product of all the irreducible factors of G taken with multiplicity 1. Also, 
D W (G) denotes the discriminant of G with respect to the variable w, i.e. the 
resultant of G and its partial derivative with respect to w. We write this 
resultant as D W (G) = Res^(G, |^). Furthermore, the following polynomials 
are introduced: 



Result: the main result of [T] is the following (see Theorem 4 and Theorem 13 



Theorem 2 Let S be an algebraic surface implicitly defined by F e R[x, y,z], 
fulfilling the above hypotheses. Then the following statements hold: 

(1) If R is not identically zero, then the set of real roots of R is a critical set 
of S. If R has no real roots, then the elements of the family show just one 
topology type. 



A. 




if deg x (M) = 

D x (M(x, z)) otherwise 



in P): 
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(2) If R is identically zero, then M = M(z), and the set of real roots of M is a 
critical set of S. 



2.2 Geometrical interpretation of the results for the implicit case. 



Let us provide a geometrical interpretation of the process, according to The- 
orem |2l giving rise to a critical set of S. For this purpose, let C = V(F, F y ) 
be the algebraic variety defined by F and its partial derivative F y (i.e. the 
intersection of the surfaces defined by F and F y , respectively). Under the as- 
sumptions made on S, one may check that gcd(F, F y ) = 1; so, C has dimension 
1, i.e. it is a space algebraic curve. Taking into account that the normal di- 
rection to the surface S at each point P e S is defined by the gradient vector 
VF(P) = (F X (P) , F y (P) , F Z (P)) , one may see that C consists of the points 
of S where the normal vector is either parallel to the xz-plane, or identically 
zero (i.e. singularities of S). Also, let us denote the curve defined by the poly- 
nomial M(x, z) on the xz-plane, as M,. Now, from a geometric point of view, 
in order to compute a critical set by means of Theorem [21 one has to perform 
two different phases: 

(1) Computation of M, (Projection Phase): from basic properties of resultants 
(see for example p. 255 in [18]), one may see that the curve M. consists of 
the projection onto the rcz-plane of C, together with the curve defined by the 
leading coefficient w.r.t. y of F, denoted as lcoeff^(F); since by hypothesis 
lcoeS y (F) does not depend on x, this last curve consists of finitely many 
lines z — z = where z is a root of lcoeff y (F). 

(2) Computation of the critical set (Analysis of the Projection): the following real 
^-values must be computed: 

(i) the z-coordinates of the points of M. with tangent parallel to the x-axis 
(included the values z a 's so that z — z a is a component of M). 

(ii) the z-coordinates of the singularities of Ai 

(iii) the values ZbS so that z — Zb is a horizontal asymptote of Ai 

In the case M = M(z) this is clear. In other case, one must notice that 
from well-known properties of resultants (again, p. 255 in [IB]), the roots 
of Res x (M,M X ) correspond either to the z-coordinates of the solutions of 
the polynomial system {M(x,z) = 0,M x (x,z) = 0}, or to the roots of 
lcoeffa;(M). Then, it suffices to interpret the solutions of the system, and 
the real roots of lcoeff x (M), from a geometric point of view. 

In the rest of the paper, we will refer to these z-values as (l)-values, (2)-values 
and (3)-values, respectively. 
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2.3 Notation and Hypotheses for the parametric case 

In the rest of the paper, let T be a family of rational algebraic curves, alge- 
braically depending on a parameter A, defined by the parametric equations 

x = u(t, A) 

< 

k y = v(t, A) 

where u, v are real, rational functions of the variables t, A in reduced form 
(i.e. the numerator and denominator of u, v share no common factor), not 
identically 0. Thus, for almost all real values of A the above equations define a 
rational curve of parameter t. In fact, the only exceptions are the values of A 
causing that some denominator of u, v vanishes. Now our purpose is to study 
the topology types arising in T as A moves in R. In our approach, the key for 
computing these topology types is the computation of a finite partition of R 
so that for each element of the partition, the topology type of the family stays 
invariant. Hence, in the sequel we focus on this question. For this purpose, 
observe that the equations 

x = u(t, A) 
< y = v(t,\) 
z = A 

define a rational surface S in parametric form, whose level curves are exactly 
the members of the family T . Hence, our aim is the computation of a critical 
set of S. For this purpose, we introduce the following notation: 

• 4>\(t) = (u(t, \),v(t, A)) is the parametrization of the family T . Hence, the 
associated surface S is defined by the parametrization (u(t, \),v(t, A), A). 

• F x (x,y) = F(x,y, A) defines the implicit equation of T . Hence, the implicit 
equation of S is F(x, y, z) = 0. 

Now in order to make precise the hypotheses that we require on the family 
J 7 , we need to recall the notion of proper parametrization. One says that a 
parametrization (p(t) of a rational curve 7i is proper, if there are just finitely 
many points of 7i generated simultaneously by several different values of the 
parameter t; intuitively speaking, this means that the curve is traced just once 
as t moves in E. The interested reader may find more information on this topic 
in Chapter 4.2 of [IS]. Then, we consider the following hypotheses: 
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Hypotheses: 



(Hi) The parametrization 4>\(t) = (u(t, A), v(t, A)) is proper for almost all values 
of A (i.e. there are just finitely many values of A such that 4>\{t) is not 
proper). 

(H 2 ) deg^Fx) = deg(Fx), i.e. considering F\ as a polynomial in the variables 
x, y, the degree of F x w.r.t. the variable y equals the total degree of F x ; 
in particular, the leading coefficient of Fx w.r.t. y does not depend on x, 
though it may depend on A. 

(H3) The function u(t, A) depends on the variable t. 

In the following, we will refer to these hypotheses as (Hi) and (H 2 ), respec- 
tively. Observe that (H 2 ) is slightly more restrictive than the hypothesis (c) 
required in the implicit case. On the other hand, (H 3 ) guarantees that Fx 
(and therefore F) depends on the variable y. Notice that if (-£^3) does not 
hold, the analysis is reduced to the family x = u(X) (consisting just of finitely 
many lines normal to the x-axis), and the problem is trivial. Finally, since S 
is a rational surface then it is irreducible; in particular, S cannot represent a 
plane normal to the z-axis because from the parametrization it is clear that 
the z-coordinate cannot be constant. Thus, whenever the above hypotheses 
hold, the surface S described by the considered parametrization satisfies all 
the hypotheses required in the implicit version of the problem. 

The problem of checking the above hypotheses (Hi) and (H 2 ) is addressed in 
the next section. Checking (H%) is trivial; so, in the sequel we assume that 
(H 3 ) holds. 



3 Checking Hypotheses 

In this section we provide an algorithm for checking the hypotheses (Hi) and 
(H 2 ) introduced in Subsection 12.31 In addition, the algorithm addressed here 
provides also a list of finitely many "special" values of the parameter A, which 
will be important (see the next section) in order to compute a critical set of 
the surface S associated with the family. 

3.1 Checking hypothesis (Hi) 

In this subsection we consider the problem of checking whether (Hi) holds, or 
not. For this purpose, we start with a technical lemma on the behavior of the 
gcd of two polynomials of R[x, y, A], under the specialization of the parameter 
A. Hence, let <p a be the natural homomorphism of R[x, y, A] into E[x, y], i.e. for 



7 



a 6 1, 



(p a : R[x, y, X] -> K[x, 

f(x,y,\) -> <p a (/) = f(x,y,a) 

Moreover, if /, # £ R[z, y, A], we write / = / • gcd(/, g),g = g- gcd(/, and 
we define the following sets: 



BP v) (f,g) = {«£ K|^a (lcoeff x (/)) = or ^ (lcoeff^)) = 0} 
Bfr v \f,g) = {ae R\<p a (lcoeff,(/)) = or ^ a (lcoeff y (#) ) = 0} 
Bi x>v \f,g) = {ae R\ Va (Res x {f,g))} = 
B%' v \f,g) = {ae R\ip a (Res y (f,g))} = 
i3(/,( 7 )=S 1 US 2 US3US 4 



Then the following lemma holds. 

Lemma 3 Let f,g £ R[x,y, A]. Then, for alia ^ B(f,g), it holds that gcd(<p a (f),<p a (g)) = 

^a(gC<l(/,y)). 

Proof. Writing f = h ■ h, g = g ■ h, it holds that 

gcd(Vo(/), ^(fiO) = gcd(^ a (/), ¥>aGj)) ' VaCO 

So, we have to prove that for a ^ g), gcd(<^ a (/), <p a (g)) = 1- Indeed, if this 
equality does not hold then either Res x ((p a (f), <p a (g)) = or Res y (<p a (f), <f a (g)) = 
0. On the other hand, since a (fi B(f, g) then the resultants behave well under 
specializations (see Lemma 4.3.1 in [21]), and hence Res x (<p a (f),<p a (g)) = 
(resp. Resy((p a (f),(p a (g)) = 0) iff (p a (Res x (f,g)) = (resp. (p a (Res y (f,g)) = 
0). Nevertheless, (p a (Res x (f , g)) = (resp. <p a (Res y (f, g)) = 0) cannot happen 
because a £ B(f,g). □ 

Now we fix the following notation: 



</> x (t) = (u{t,X),v(t,X)) 



( 



X n (t,X) X 21 (t,X) 

x 12 (t,xy x 22 (t,x) 



) 



and we introduce the following polynomials: 



Gf(t, s, A) = X u (t, A) • X 12 {s, A) - X 12 (t, A) ■ X 11 (s, A) 
Gl(t, s, A) = X 21 (t, A) • X 22 (s, A) - X 22 (t, A) • X 21 (s, A) 
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Moreover, we write G\ = G\ • gcd(Gf,Gf), G% = Gf ■ gcd(Gf,Gf), and we 
denote: 

V\ = B?' s) (Gf, Gt), V 2 = B i 2 t ' s) (Gf, Gt) 
V 3 = Bt\Gt,Gl), V, = Bf s) {GiGt) 
D = DiUD 2 UD 3 U £> 4 , 

(notice that T>\ = V 2 , T> 3 = U± because Gf,Gf are symmetric w.r.t. t,s). 
Then, the following theorem holds. Here, we denote the evaluations of Gf, Gf 
at A = a, as ip a (Gf), ip a (G 2 ), respectively. 

Theorem 4 The parametrization 4>\(t) is proper for almost all values of X iff' 
deg f (gcd(Gf, Gf)) = 1, i.e. iff gcd(Gf , G 2 ) =t — s (maybe multiplied by some 
polynomial a(\) G K[A]j. Furthermore, if this condition holds, then the only 
values of X where 4>\(t) may not be proper, are those in T>. 

Proof. Let a£i By Theorem 4.30 in [18], <f) a (t) is proper iff 
deg, (gcd (MGt),MGt))) =1. 

However, by Lemma El if a ^ T> then 

gcd (MGf), <p a (Gti) = Va (gcd (Gf, Gt)) 

Hence, if we prove that 

deg t (<p a (gcd (Gf , Gt)))= deg t (gcd(Gf , Gf)) 

for a £ T>, then we have finished. Indeed, if this equality does not hold, then 
the leading coefficient w.r.t. t of gcd(Gf , Gf) vanishes at A = a. However, in 
that case the leading coefficients of Gf, Gf w.r.t. t both vanish at A = a, and 
this cannot happen because a £ T>. Finally, notice that t — s always divides 
gcd(Gf,Gf); hence, deg t (gcd(G"f, Gf)) = 1 iff gcd(Gf,Gf) =t-s (perhaps 
multiplied by some polynomial a (A)). □ 

Hence, Theorem H] gives us an algorithm for checking hypothesis (Hi). More- 
over, if the condition in Theorem H] holds one can determine the set T>, which 
contains the finitely many values of the parameter where properness fails. If 
the condition does not hold, one can compute a reparametrization £\(t) of the 
family, proper for almost all values of A, by applying the reparametrization 
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algorithm in Section 6.1.2, p. 193 of [18]. We give more details in Subsection 

El 



Remark 1 Observe that, if they exist, the X-values making that the denomi- 
nator of either u or v is identically belong to D . 

3.2 Checking hypothesis (H 2 ) 

Now let us consider hypothesis (H2). For this purpose, we recall that the 
degree of a rational function (i.e. of a quotient of polynomials) is defined as the 
maximum of the degrees of the numerator and the denominator; furthermore, 
the degree of a rational parametrization 4>(t) = (xi(t), is defined as the 

maximum of the degrees of Xi{t)i X2{t) (which are rational functions). We also 
recall the following result from [18] (see Theorem 4.21 in |18j). 

Proposition 5 If <p(t) = {xi{t) , X2{t)) , where Xi{t) i> s n °t identically 0, is 
a proper rational parametrization of a curve and f(x,y) = is its implicit 
equation, then deg t (xi(t)) = deg y (/) 

Moreover, we also need the following lemma (see for example Section 3 of 

to. 

Lemma 6 Let f(x,y) implicitly define a plane curve V, let \x G R, and let 
9^{ x i y) be the implicit equation of the curve that is obtained from V by applying 
the linear transformation {x = X + (jY,y = Y}. Then, for almost all values 
of n it holds that deg y (g M ) = deg(g M ) (here, deg^) denotes the total degree 

ofgj- 

Now the following result holds. This theorem provides a method for computing 
the degree of a rational curve just from its parametrization, without making 
use of its implicit equation (compare also with Theorem 6 and Theorem 7 in 

m- 

Theorem 7 Let <p(t) = (xi(t), be a proper rational parametrization of 

a curve V, where Xiif) ^ s n °t identically 0; also let Xn(t) — Xiif) ~ f i X2{t), 
with fj, generic. Then, deg(V) = deg t (x M (t)). 

Proof. Let be the curve obtained from V by applying the linear trans- 
formation {x = X + fjY,y = Y}. Then, for all values of \x it holds that 
is a rational curve properly parametrized by M (t) = (Xn(t)jX2(t))- Now let 
g^{x,y) be the implicit equation of Up. Then by Lemma for a generic fi 
it holds that deg y (g^) = deg(g fl ). Now since the degree of a curve is invari- 
ant by linear transformations, we have that deg(V) = deg^). Finally, since 
Xi(£), X2(t) are not both identically 0, for a generic /1 it holds that Xn(t) * s n °t 
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identically 0, either; moreover, 4>n(t) is proper for every \i G R, and hence by 
Proposition [5] we have that deg y (g^) = deg t (% /i (t)). Therefore the statement 
follows. □ 

Theorem [7] together with Proposition [5] provides the following corollary. 

Corollary 8 Let <f)(t) = (xi(t), X2(t)), where Xi(t) n °t identically 0, be a 
proper rational parametrization of a curve V , and let f(x, y) = be its implicit 
equation. Also, letx^it) = Xi(t) ~ l l X2(t) , with [i generic. Then, deg t (xi(t))) = 
deg t (x M (t)) iffdeg y (f) = deg(f). 

Then the following characterization of {H 2 ) can be deduced. 

Corollary 9 {H 2 ) is fulfilled iff, for a generic [i, it holds that deg t (u(t, A)) = 
deg t (u(t, A) — fiv(t, A)). Moreover, the above value provides the degree of F\, 
as a polynomial in the variables x,y. 

Also from [10], one may see that if (-f^) is not satisfied, almost all changes 
of coordinates of the type {x = X + /iY, y = Y} set the surface properly. 
Furthermore, if (-£^2) is fulfilled, then the following proposition holds. 

Proposition 10 Assume that (H 2 ) holds, and let m = deg t (Xu(t, X)) , n = 
deg t (X 12 (t,X)), a(A) = \coeS t (X n (t, X)), 6(A) = \coeS t (X 12 (t, A)). Also, let 
A be a real root of the leading coefficient of F\(x,y) w.r.t. the variable y, so 
that: (i) Xn(t, Xq) does not vanish identically; (ii) <p\ {t) is proper. Then, the 
following statements are true: 

(1) Ifm>n, then a(Ao) = 0. 

(2) Ifm = n, then a(A ) = b(X ) = 0. 

(3) Ifm<n, then b(X ) = 0. 

Proof. Since by hypothesis (-^2) holds, if the leading coefficient of F\(x,y) 
vanishes at A = A then deg y (F\ ) < deg y (F\). Also, because of the re- 
quired conditions (i) and (ii), Proposition [5] holds and therefore deg y (F\ a ) = 
deg t (xi(t, A )). Moreover, from Theorem[7Jit holds that deg y (F\) = deg t (xi(£, A)). 
So, we deduce that deg t (xi(£, A )) < deg t (xi{t, A)). The rest follows from the 
definition of degree of a rational function. □ 

Notice that if Xu(t, Xq) = 0, then Ao is an element of the set D in Theorem 
HI Now Proposition [10] provides the following corollary. Here, we keep the 
notation used in the above proposition. 

Corollary 11 Assume that {Hi) and (H 2 ) are fulfilled. Then, the real roots 
of\coeS y (F) belong to the set consisting of the following elements: (i) the real 
elements of the set V in Theorem^ (ii) the real roots of: (a) a(X), if m > n; 
(b) b(X), if m < n; (c) gcd(a(A), b(X)), if m = n. 
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3. 3 Normality 



Given a parametrization i/j(t) of a plane curve S, one says that the parametriza- 
tion is normal if it is surjective, i.e. if for all Pq G £ there exists ^ G C so that 
if) (to) = Po] notice that to m ay be complex even though P Q is real. This notion 
has been studied in [2], [17], [18]. For our purposes, here we recall the following 
result (see Theorem 6.22 in |18j). 

Theorem 12 Let (bit) = ( — 11 . . , — | be a parametrization of a plane 

vyi \x 12 (ty x 22 (t) J 

curve £, and let n = deg(X 12 ), s = deg(X 22 ). Also, let b* = coeff (Xu, ra) (i.e. 
the coefficient of degree n in Xn(t)), b = coeff (X 12 , n), d* = coeff (X 2 i, s), 
d = coeS(X 22 , s). Then, 

(i) If m > n or r > s, then <p(t) is normal. 

(ii) If m <n and r < s, then (p(t) is normal iff 

de gi (gcd (b*X 12 (t) - bX u (t), d*X 22 (t) - dX 21 (t))) > 1 

Moreover, if <p(t) is not normal, the only point that is not reached by the 

(b* d*\ 

parametrization is the so-called "critical point" — , — ; which is a point 



b ' d 



ofS. 

In our case, we deal with the family 

'x n (t,\) x 21 (t,xy 



<t> x (t) = (u(t,X),v(t,X)) 



x 12 (t,xy x 22 (t,\) t 



Hence, from Theorem [T2l we can derive the following result: 

Theorem 13 Letm = deg f (X n ), n = deg t (X 12 ), r = deg t (X 2 i) ; s = deg t (X 22 ). 
Also, let a(X) = coeS(X n (t, X),m), b*(X) = coeff (X n (t, X),n), b(X) = coeS(X 12 (t, X),n), 
c(X) = coeS(X 21 (t,X),r), d*(X) = coeS(X 21 (t,X),s), d(X) = coeff (X 22 (t, A), s). 
Then, it holds that: 

(i) If m > n and r > s, the values of X where <fi\(t) may not be normal satisfy 
gcd(a(A),c(A)) = 0. 

(ii) If m > n and r < s, the values of X where <j)\(t) may not be normal satisfy 
a(X) = 0. 

(Hi) If r > s and m < s, the values of X where 4>\(t) may not be normal satisfy 

c(A) = 0. 
(iv) If m < n, r < s, and 

deg, (gcd(b*(X)X 12 (t)-b(X)X 11 (t),d i '(X)X 22 (t)-d(X)X 21 (t))) > 1, 
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the values of X where <p\(t) may not be normal satisfy at least one of the fol- 
lowing conditions (here, we denote rj(t, A) = \coeS t (b*Xi 2 — bXu), u(t, A) = 
d*X 22 - dX 21 , fj = 7}/ gcd(?7, v), v = vj gcd(?7, v)): (a) lcoeS t (r](t, A)) = 0; 
(b) lcoeS t (u(t, A)) = 0; (c) Res t (fj, v) = 0. 
(v) If n < m, r < s, the only points of the surface S that may not be reached 
by the parametrization, are the points of the space curve C cr u parametrized 

, b'(\) d*(\) \ 

Proof. The statements (i), (ii), (iii) and (v) essentially follow from Theorem 
[l"2l Statement (iv) follows from Theorem [T2l and Lemma 4.26 in [18]. □ 

3.4 Summary, and an example 

Here we provide the full algorithm for checking the hypotheses (Hi), (H 2 ), re- 
quired on the family T , and we illustrate it by means of an example. Besides 
checking (Hi) and (H 2 ), the algorithm also provides a list of finitely many 
'special" values of the parameter, which will be useful in the next section. 
If (Hi) does not hold, then one can use the reparametrization algorithm in 
Section 6.1.2, p. 193 of [18] in the following way (see [18] for further infor- 
mation): (1) let = gcd(Gf, Gt); then, choose (a, (5) G K 2 so that G 4 '(a,(5) 
is not identically (those finitely many A-values making G^(a,f3) = are 
incorporated to the list of "special" values); choose also a,b,c,d G M so that 
ad — be 7^ 0. (2) Consider the rational function 

= aG*(a,t) + bG<l>(M 
AU cG*(a,t) + dG*(P,t) 

Then, let r = deg((fi\) / deg(R\(t)) . In general the value of r has to be dis- 
cussed upon the value of A, but it will be constant except for finitely many 
A-values, which again must be stored in the list of special values. (3) Introduce 
a generic rational parametrization Q(t) of degree r (i.e. the generic value of r) 
with undetermined coefficients. From the equality <fi\(t) = Q(R\(t)) derive a 
linear system of equations in the undetermined coefficients, and by solving it 
determine Q(t) (in fact Q\(t)). Notice that this linear system has A as a pa- 
rameter, and therefore, again, certain (special) values of the parameter must 
also be computed (by discussing the system). In the end, the whole process 
yields a reparametrization (depending on A), and a finite list of special values 
of A. The correctness of this process follows from Theorem 6.4, p. 191 of [18]. 

Also, if (-^2) is not satisfied almost all changes of coordinates of the type 
{X = x + fiy, Y = y} set the surface properly; observe that this kind of 
transformation leaves the z-coordinate invariant, and so the topology of the 
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level curves of the surface is not changed. 

Algorithm: (Check) Given a uniparametric family T of rational curves, defined 
by its parametric equations 



where m = deg t (X n ), n = deg t (Xi 2 ), r = deg t (X 2i ), s = deg t (X 22 ), and 
a(X) = coeff(X n (t,A),m),&*(A) = coeff (X n (t, X),n),b(X) = coeft(X 12 (t, A), n), 
c(A) = coe&(X 2 i(t,X),r),d*(X) = coeff (X 21 (t, A), s), d(A) = coeff (X 22 (t, A), s), 
and where A is a real parameter, the algorithm checks hypotheses (Hi) and 
(H 2 ); moreover, if the hypotheses are satisfied, the algorithm also computes a 
list Spec of finitely many "special" values of A. 

1. (Check hypothesis (Hi)) 

1.1 Compute the polynomials Gf,G 2 in Section I3TT1 

1.2 If deg t (gcd(Gi , G 2 )) = 1, then return (Hi) holds, otherwise return (Hi) 
does not hold. 

1.3 If (Hi) does not hold, reparametrize the family. Let Spec be the list of 
special values computed in the process (if (Hi) holds, Spec = 0). 

1.4 Compute the set T> in Theorem HI and let Spec x := T>. 

2. (Check hypothesis (H 2 )) 

2.1 Compute deg t (u(t, A)) , deg t (u(t, A) — ^v(t, A)). If both are equal, then 
return (H 2 ) holds, otherwise return (H 2 ) does not hold. 

2.2 If (H 2 ) does not hold, apply a change of coordinates {X = x + fiy, Y = y}, 
and go back to 1.3. 

2.3 If (H 2 ) holds, then compute the real roots of: (a) a(X), if m > n; (b) b(X), 
if m < n\ (c) gcd(a(A), b(X)), if m = n. Let Spec 2 be the set consisting of 
these values. 

3. (Normality) 

3.1 If m > n and r > s, let Spec 3 be the set of real roots of gcd(a(A), c(A)). 

3.2 If m > n and r < s, let Spec 3 be the set of real roots of a(X). 

3.3 If r > s and m < n, let Spec 3 be the set of real roots of c(A). 

3.4 If If m < n, r < s, then compute 

8 = deg t (gcd (b*(X)Xi 2 (t) - b(X)Xn(t), d*(X)X 22 (t) - d(X)X 2 i(t))) 

3.4.1 If 5 > 1, then let r/(t, A) := lcoeff 4 (6*X 12 - 6X n ), v(t, X) := d*X 22 - 
dX 2 i, f] := 77/ gcd(r/, u), v := z// gcd(r/, z/), and let Spec 3 be the set of 
real roots of lcoefft(?7(£, A)), lcoefft(^(t, A)), Res t (fi, v) = 0. 

3.4.2 If 5 < 1 then determine the parametrization 




b*(X) d*(X) 



A 



6(A) ' d(X) 
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of the space curve C cr u. Let Spec 3 : = 
4. Special values: Spec := Spec Q U Spec ± U Spec 2 U Spec 3 



The above algorithm is illustrated by the next example. Here, we consider 
the parametric equation of the offset family to a cardioid. By using the re- 
sults in j3], one may check that this offset family is rational; moreover, the 
parametrization used in the example is taken from [T6] . 

Example 1 Let Od{C) be the offset family to the cardioid (x 2 + Ay + y 2 ) 2 — 
16(x 2 +y 2 ) = 0. Then, a parametrization of this family is <fid{t) = (u(t,d),v(t,d)), 
where: 



Here, the parameter d denotes the offsetting distance. So, let us apply the 
algorithm Check. In step (1), we analyze (Hi). For this purpose, we compute 
Gf,Gt (step 1.1): 

Gf = (3456t 5 -31104t 3 +rft 8 -126rft 6 +10206rft 2 -6561d)(486s 4 +36s 6 +2916s 2 + 
s 8 + 6561) - (3456s 5 - 31104s 3 + ds 8 - 12Qds 6 + 10206ds 2 - 6561d)(486t 4 + 
36£ 6 + 2916t 2 + t 8 + 6561) 



G% = -18t (864t 3 - 16t 5 - 1296t + dt 6 - 21dt 4 - 189dt 2 + 729d) (486s 4 + 36s 6 + 
2916s 2 + s 8 + 6561) + 18s(864s 3 - 16s 5 - 1296s + ds 6 - 21ds 4 - 189ds 2 + 
729d)(486t 4 + 36t 6 + 2916t 2 + t 8 + 6561) 

Then, we check that gcd(Gf , Gf) = t—s, and therefore that deg t (gcd(Gf , G 2 )) = 
1 (step (1.2)); hence, (ifi) holds, Spec := 0, and we go to step (1-4)- Here, 
we compute the set T> containing the d-values where <frd{t) may not be proper; 
for this purpose, we compute the real roots of: (i) Content s (\coeS(Gf,t)); (ii) 



Content s (lcoeS (G 2 , t)); (Hi) Content^ (Rest(Gf , G 2 ) ) . In this case, we get that 



Content s (\coeS(Gf,t)) = 54, Contents (lcoeff(Gt t)) = 18, andRes t (Gf, G 2 ) = 
C ■ d{t 2 + 9) 44 (729d - 1053dt 2 - 117dt 4 + dt 6 + 3456t 3 ) (-32* + dt 2 + 9d) 2 (with 
C e Nj. Hence, it holds that D = {0} and therefore 

Spec 1 := {0}. 

Now in step (2), we check (H 2 ). For this purpose, in step (2.1) we compute 
deg t (u(t,d)) = 8 and deg t (u(t,d) — fiv(t,d) = 8; since both of them coincide, 
then (H 2 ) holds, and we go to step (2.3). Since m — n — 8, we compute 



u(t, d) 



v(t, d) 



3456t 5 - 31104t 3 + dt 8 - 126dt 6 + 10206dt 2 - 6561rf 

486i 4 + 36t 6 + 2916t 2 + t 8 + 6561 
(-18)t(864t 3 - 16t 5 - 1296* + dt 6 - 21dt 4 - 189dt 2 + 729d) 

486t 4 + 36t 6 + 2916* 2 + t 8 + 6561 
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a(d) = d, b(d) = 1, and we determine gcd(a(d) , b(d)) , which is 1; hence, 
S r pec 2 := 0. 

Finally, in step (3) we consider normality questions. We have thatr = 7, s = 8 
and therefore r < s. Since m = n, we go to (3.4) and we compute 

7}{t, d) = c/(486t 4 + 36t 6 + 2916t 2 + t 8 + 6561) - 3456t 5 + 31104t 3 - dt 8 + 126dt 6 

-10206dt 2 + 6561c/ 
u(t, d) = 18c/t(864t 3 - 16t 5 - 1296t + dt 6 - 2ldt 4 - 189dt 2 + 729d) 
^ = deg t (gcd(?7, v)) 

We get 5 = 0; hence, we go to step (3.4-2) and we obtain a space curve of 
possibly non-reached points, namely 

Ccrit = (d, 0, d) 
Finally, we get that Spec := {0}. 



4 The algorithm for the parametric case. 

Based on the observations made in Subsection 12.21 in this section we present 
an algorithm for computing a critical set of S, under the assumption that 
the hypotheses requested in Subsection 12.31 are satisfied; the algorithm does 
not compute or make use of the implicit equation of F. More precisely, in 
Subsection 12.21 it is observed that one can compute a critical set of S by 
determining certain z-values referred to as (l)-values, (2)-values and (3)-values, 
respectively. These values have different geometric meanings, and are related 
to notable points and lines of the curve M. defined in the icz-plane by the 
polynomial M(x, z) = \jD y (F), where F is the implicit equation of S. On the 
other hand, also in Subsection 12.21 it is shown that A4 can be written as the 
union of the projection onto the rrz-plane of the space curve C = V(F,F y ), 
and the curve (in the xz-plane) defined by \coe& y (F). Now, the real roots of 
lcoeS y (F) (or, more precisely, a finite set containing them) can be determined 
by means of Corollary [Til without explicitly computing F. Hence, in this 
section we focus on the computation of the remaining (1), (2) and (3)-values, 
which are related to C; for this purpose, and since we do not want to compute 
or make use of F, the idea is to work with a "parametric description" of C. 
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In order to provide equations for C, one may see (recall the first paragraph in 
Subsection 12. 2j) that this curve consists of the following points: 



• points of S, reached by the parametrization, such that the normal vector to 
S is either (in which case the point is a singularity of S), or parallel to 
the xz plane. 

• self-intersections of S, reached by the parametrization. 

• points of S with some of the above geometric properties, but not reached 
by the parametrization. 

In the sequel we will refer to these sets as first, second and third sets, respec- 
tively. Notice that some of these sets may be empty, and that they are not 
necessarily disjoint. So, we start with the first one. By performing easy com- 
putations with the parametrization (u(t, X),v(t, A), A) of S, one may see that 
the expression of the normal vector to S is: 

N = v t i — u t j + (u t v x - u x v t )k 



Hence, we define C\ as the following subset of points (x, y, z) e 



-3. 



Ci 



x = u(t, A) 
V = v(t,\) 
z = A 
h{t,X) = 

X 12 (t,z)-X 22 (t,z) ± 



where h(t,X) is the square-free part of the numerator of u t (t,X). One may 
see that C\ contains the first set. Then let us consider the second set (i.e. 
self-intersections of S). For this purpose, we impose 

u(t, A) = u(s, X),v(t, A) = v(s, A), with t ^ s 



Hence, we find again the polynomials Gf, G 2 introduced in Subsection 13.11 
We recall from Subsection 13.11 the notation Of , G 2 for the polynomials ob- 
tained by removing the common factor t — s in Gf, G 2 . Furthermore, we write 

j(t, A) = \jRes s (Gf , G 2 ). Since by hypothesis the parametrization 4>\(t) = 
(u(t, A), v(t, A), A) is proper for almost all values of A, then j(t,X) cannot be 
identically (see Theorem 4.30 in [IS]) - Therefore, we define C 2 as the following 
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subset of points (x, y, z) G 



Co 



x = u(t, A) 

y = v(t,\) 

z = A 
j(t,\)=0 

X 12 (t, z)-X 22 (t, Z)yt0 



One may see that this set contains the self-intersections of S reached by the 
parametrization. Finally, the third set of points has been studied in Subsection 
13.31 In this sense, by Theorem [TBI we have basically two possibilities: in the first 
case (see statements (i), (ii), (iii), (iv) in Theorem [T31 . there exists a finite set 
TV containing the A- values so that 4>\(t) may not be normal; in that situation, 
we just add TV to the rest of elements of the critical set, and therefore no 
special difficulty arises. In the second case (statement (v) of Theorem [T3|) . 
it may happen that there are infinitely many A values so that <p\ (t) is not 
normal; thus, the non-reached points of 5* are among the points of the space 

'6*(A) d*(\) 



curve C cr u parametrized by 



b(\) ' d(X) 

assume that we are in this more complicated situation. 



A . Hence, in the sequel we will 



4-1 Computation of a critical set 



Since the roots of lcoeS y (F) can de determined from Corollary [TT], we focus 
on the remaining elements of the critical set. For this purpose, we need to 
analyze the projection of C onto the xz-p\&ne. From the above reasonings, we 
can write 



C C Ci U C 2 U C 



crit 



So, n xz (C) C n xz {Ci) U 7T XZ {C 2 ) U 7T xz {C crit ). Now we have that 



x = u(t, z) 
h(t, z) = 

X 12 (t,z)-X 22 (t,z)^0 



Let f(x, t, z) be the numerator of x — u(t, z). Also, we write h = h-gcd(h, X\ 2 ■ 



X 22 ). Then, we denote the curve defined by m^(x, z) = \fRes t (f, h) as M. 



Hi 
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notice that m^(x, z) cannot be identically because by hypothesis the func- 
tion u(t, A) is given in reduced form (i.e. its numerator and denominator share 
no common factor). Furthermore, one may see that tt xz (Ci) is included in A4\. 
Similarly, 



n xz {C 2 ) = < 



x = u(t, z) 
j(t,z) = 

X 12 (t,z)-X 22 (t,z)^0 



In this case, we write j = j ■ gcd(h,Xi 2 ■ X 22 ), we denote m^ 2 '(x, z) = 
Kes t (f,j) (which cannot be identically for the same reason than m^(x, z)), 



and we represent the curve defined by rrS 2 ^ as Ai 2 . One may observe that 
^xz(C 2 ) is included in A4 2 . Finally, 7r X2 (C cr j t ) is the parametric curve defined 

Now, the ^-coordinates of the singularities and of the points with tangent 
parallel to the re- axis of M.\ (resp. A4 2 ), and also the values z a so that z — 
z a is a horizontal asymptote of M.i (resp. -M 2 ), correspond to real roots of 
Res x (m^), m^) (resp. Res x (m^ 2 \ m^)). On the other hand, since 7i xz (C crit ) 
f b* (X) \ 

is parametrized by ( . , Aj , it has no point with tangent parallel to the 

x-axis, no singularity, and no component parallel to the x-axis. However, it 
may have horizontal asymptotes, corresponding to the roots of b(X). 

Finally, in order to compute a critical set we also need to determine the inter- 
sections between M. 1 and M. 2 , M.i and 7i xz (C crit ), M. 2 and 7i xz (C crit ), respec- 
tively. In the first case, let 

777,(1-) _ m (i) j gcd(m ( - 1 ' ) , m ( - 2 ' ) ), m( 2 ) = rrS 2 ^ j gcd(m^\ rrS 2 ^); 

then, the ^-coordinates of the intersections between A4i and A4 2 are contained 
in the set of roots of Res x (m( 1 \ m^). In the second (resp. the third) case, we 
simply compute the roots of the numerator M^'(z) (resp. M^^z)) where 

M^(z) (resp. M^ 2 \z)) is defined as the result of substituting x := , 

6(A) 

z := A in m^(x,z) (resp. m^(x, z)), whenever this substitution does not 
yield 0, and is defined as 1 in other case. 

So, we can derive the following algorithm for computing a critical set of T . 
Here, we will use the subset Spec computed by the algorithm Check, that we 
developed in Section [3J 
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Algorithm: (Critical) Given a uniparametric family T of rational curves de- 
pending on a parameter A, defined by its parametric equations 

fulfilling hypotheses (#1), (H 2 ), (H s ), where m = deg t (Xn), n = deg t (Xi 2 ), 

r = deg t (X 2 i), s = deg t (X 22 ), and 6*(A) = coeS(X n (t, A), n), 6(A) = coeff (X 12 (t, A), n), 

d*(X) = coeff(X 21 (t, A), s), d(X) = coeS(X 22 (t, X), s), the algorithm computes 

a critical set A of the family. 

(1) Compute the set Ai consisting of the real roots of the following polynomials: 

• Res^m*- 1 ), mW) 

■ Res x (mW, m^). 

(2) If m > n, or r > s, or m < n, r < s but 

deg t (gcd(6*(A)X 12 (t) -b(X)X 11 {t),c{X)X 22 {t) - d{X)X 21 {t))) > 1, 

then A 2 = 0. Otherwise, let A 2 be the set consisting of the real roots of the 
following polynomials: 

■ 6(A). 

• AfW(z) 

• M( 2 )(z) 

(3) Let A = Spec UAiUA 2 . Return A. 
4-2 Correctness of the Algorithm 

The aim of this subsection is to prove that the algorithm Critical, provided 
in the above subsection, is correct, i.e. that the set A determined by the 
algorithm is a critical set of S. The necessity of this proof is due to the fact 
that in Critical we are working not really with ir xz (C), but with the curve 
Jli* = JH 1 \j J\4 2 U n xz (C cr it), which in general may be bigger than n xz (C). 
Hence, we want to ensure that no critical value (i.e. no A-value where the 
topology of the family changes) has been missed when passing from tt xz (C) to 
At*. This is done in the following theorem. 

Theorem 14 The algorithm Critical is correct. 

Proof. Since tt xz (C) C M* , from Chapter 3, Section 2.5, Exercise 5 of [2U] it 
follows that the singularities of ir xz (C) are also singularities of Ai*, and hence 
their ^-coordinates are found by the algorithm. Now, let Q G n xz (C) so that 
the tangent line I to tt xz (C) at Q is parallel to the x-axis. If Q is a regular point 
of Ai*, since ir xz (C) C Ai* we have that £ is also the tangent line to Ai* at 
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Q; therefore Q is a point of A4* with tangent parallel to the rc-axis, and hence 
its 2-coordinate is found by the algorithm. Otherwise, Q is a singular point 
of Ai* and therefore its ^-coordinate is also determined by the algorithm. On 
the other hand, asymptotic branches of ti xz (C) and components lying in planes 
normal to the z-axis, keep their nature when passing to A4*. Finally, since the 
real roots of \coe& y (F) are included in Spec, we conclude that all the (1), (2) 
and (3)-values are computed by Critical, and the result follows. □ 

Remark 2 Notice that we are stating that the set computed by the algorithm 
Critical is a critical set, i.e. that it contains all the \-values where the topol- 
ogy type of the family may change. But we are not stating that the set that 
we are computing is optimal. So, the output of the algorithm may include 
additional z-values where the topology of the family does not change. This 
drawback was already present in the algorithm of [7]/, which was not optimal, 
either. However, our algorithm may yield critical sets which are bigger than 
those in fj]/. Since ajfine transformations preserve the topology of curves, one 
can recognize some superfluous values in a computed critical set in the follow- 
ing way: (1) apply a random afflne transformation {X = x + fiy, Y = y} to the 
family; (2) compute a critical set of the new family; (3) discard those values 
of the original critical set that do not belong to the new critical set. 



4-3 Experimentation and Results 

The algorithm Critical has been implemented in Maple, and tested with 
several examples. In this subsection, first we continue the analysis of the topol- 
ogy types in the offset family to the cardioid, started in Example [TJ Then we 
provide a simple example illustrating the non-optimality of the algorithm, as 
mentioned in Remark [2J Finally we present a table comparing timings between 
the algorithm for the implicit case, deducible from [TJ, and our algorithm. 

Example [TJ (continued): By applying the algorithm Critical, one deter- 
mines the following critical set: 

A = {-16/3, -a, -3\/3, -8^/3, -3^3/2, 0, 3^/2, 8^/3, 3\/3, a, 16/3} 

(which coincides with the output of the implicit algorithm), where a is the 
real root of 729A 5 - 1215A 4 + 702A 3 - 18A 2 + 13A - 27. The total amount of 
time required for this computation was 1.5 seconds (the cost of checking the 
hypotheses is included). From this critical set, one may deduce that there are at 
most 19 different topology types in the family. However, by applying a random 
linear transformation as suggested in Remark^, one can compute a reduced 
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critical set, namely 

{-16/3, -3V3, 0, 3^, 16/3} 



In this case, because of the properties of offset curves, one has that for d = 
one gets the original curve, i.e. a cardioid, and that for do and —do, the shape 
is the same. So, in our analysis we have just considered positive values of d. 
In Figure 2 one may find the different shapes arising in the family, and the 
intervals corresponding to each of them. In the first row (at the top of the 
figure), we display the pictures (i), (ii), (Hi) corresponding to the distances 
d = 1, d = 3, d = 22/5, respectively, all of them belonging to the interval 
(0,3\/3/2) and therefore sharing the same topology type (we have plotted the 
three pictures so as to clearly see the evolution of the family, as d is increased). 
In the second row, from left to right we have the picture (iv) corresponding to 
d = 3v^3, the picture (v) corresponding to a distance d E (3v3, 16/3), and (vi), 
that corresponds to d = 16/3. In the third row, the shape (vii) corresponding 
to d > 16/3 is shown. Also, in each figure we have included the plotting of the 
original cardioid. One may see that the offsets exhibit two cusps for d < 3y/3, 
and a loop for d > 16/3. However, the topologies of (iv) and (v) are not 
completely clear, since the picture does not show well enough the behavior 
next to the singularity with negative y-coordinate. If one enlarges the part of 
the curve next to this singularity, one obtains the pictures in Figure 3. Here 
we have plotted a detail of (iv) (left), of (v) (middle), and of (vi) (right). So, 
in (iv) there is a non-ordinary singularity; in (v) there is not one, but two 
singularities, corresponding to two self-intersections of the curve, giving rise 
to two different loops; in (vi), the topology changes so that the curve has only 
one loop (the origin, in this case, is a singular point). 

Example 2: Consider the family <fid{t) = (u(t,d),v(t,d)), where: 

u(t, d) = -25 + lit 2 - 29t + c/(57 - 95t 2 - 22t) 
v(t, d) = 49 + 18t 2 + 51* + d(70 + 34t 2 - 64t) 



The implicit equation of the family is 1136239 — 393995c? — 19629x — 53165?/ + 
2028856?x+1305309926? 3 -1200232G? 2 x+374269G??/-2090ci|/ 2 +121|/ 2 +59360320G? 4 + 
324x 2 - 396yx - 33513124c/ 2 + 1156dV + 1224dx 2 - 688992c/ 2 y - 1781936?/c/ 3 + 
9025c/ V + 2Q72dxy + QAQ0yd 2 x - 146992xc/ 3 = 0. 

Now by applying the algorithm in we get that I — j is a critical set. How- 
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ever, Critical yields ^7"} ( w ^ c ^ ^ n ^ s case coincides with Spec x , see 
Subsection 3.4}) ; so, the second element of this last set is clearly superfluous. 
In fact, when plotting curves corresponding to d G (—00, 11/25), d = 11/25 
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Fig. 2. Offsets to the cardioid 

and d £ (11/25, oo) , one gets a parabola in all the cases; hence, the topology 
type does not change for d £ R, i.e. even the value 11/25 provided by the 
implicit algorithm, is superfluous. The total amount of time required for the 
whole computation is 0.094 seconds. 

Comparison Table. The following table shows a comparison between the 
algorithm derived from the results in PQ, and our algorithm, both of them 
implemented in Maple. For each family, in this table we include whether the 
parametrization is rational or polynomial (Type: R, rational, or P, polyno- 
mial), and we provide the following data: the degree of the implicit equation 
of the associated surface (deg(/)), the total degree of the parametrization 
(deg 4 (0)) i.e. the greatest power of the parameter t arising in the numerators 
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Fig. 3. Details of Some Offsets to the cardioid 



and denominators of the coordinates, the highest power (deg A (0)) of A arising 
in the numerators and denominators of the coordinates, the timing for the 
algorithm in [T] (Imp., in seconds), the timing for our algorithm (Crit., in 
seconds), the size of the critical set determined by the parametric algorithm 
(Size), and the difference (Dif.) between the sizes of the critical sets provided 
by both algorithms (i.e. the size of the critical set provided by Critical mi- 
nus the size of the critical set provided by the algorithm in [I]). The symbol 
* in the column of Imp. means that the algorithm has been unable to provide 
an answer, or that the computation time exceeded reasonable bounds. In Ap- 
pendix I one may find the expressions of all the parametrizations used here. It 
is worth mentioning that except for the families numbers 3 and 7, the rest of 
the examples have been randomly generated. Also, the timings given include 
the cost of checking the hypotheses. 
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Family 


Type 


deg(/) 


deg t (0) 


deg A 0) 


Imp. 


Crit. 


Size 


Dif. 


1 


P. 


8 


4 


1 


35.347 


1.045 


14 


1 


2 


P. 


10 


5 


1 


160.664 


110.076 


36 


1 


3 


R. 


22 


7 


4 


* 


252.390 


70 


* 


4 


R. 


9 


6 


1 


* 


35.990 


33 


* 


5 


P. 


10 


5 


1 


158.140 


108.389 


36 


1 


6 


P. 


28 


5 


4 


39.983 


14.945 


21 





7 


R. 


16 


5 


2 


* 


31.543 


33 


* 


8 


R. 


9 


3 


2 


9.157 


0.344 


12 


3 


9 


R. 


24 


6 


3 


* 


74.770 


78 


* 


10 


R. 


6 


2 


2 


0.109 


0.093 


7 


2 


11 


R. 


12 


3 


3 


19.580 


0.389 
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12 


R. 


18 


3 


3 


66.253 


1.325 


30 


9 



4-4 Improvements in the computation. 

The (l)-values, (3)-values and (2)-values not corresponding to self-intersections 
of M. can be more efficiently computed by taking advantage of certain geo- 
metric properties of C. In fact, one can determine these values by solving 
polynomial systems in two variables; so, we can avoid one resultant, and iden- 
tify quite fast some values as potentially critical. This is based on two classical 
results. The first one follows essentially from Proposition 3 of [7]. 

Proposition 15 Let Q = (x q ,z q ) be a singularity of Ai, which is not a self- 
intersection of M., and such that z q is not a root of\coeS y (F). Then, one of 
these two possibilities occur: (i) Q is the projection of a singularity of C; (ii) 
there exists a point of C, projecting onto Q, so that the tangent to C at this 
point is normal to the xz-plane. 

The second result relates the non-singular points of At with tangent parallel to 
the z-axis, to certain notable points of C. It can be easily proven by reasoning 
with places. 

Proposition 16 Let Q G Ai be a non-singular point of Ai with tangent 
parallel to the x-axis. Then, there exists some point Q' e C, projecting onto 
Q, so that the tangent to C at Q' is parallel to the xz-plane. 
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So, let us consider first (1)- values and (2)-values. Those of these values not 
corresponding to: (i) self- intersections of Ai, (ii) real roots of lcoeff y (F), 
(hi) points of n crit (C), are real z -va\ues fuhlling that there exists (x , z ) e 
A4i U M.2- Now Mi can be seen as the union of the following two curves: (1) 
the projection onto the xz-plane of the space curve C\ defined by f(x,t, z) = 
0, h(t, z) = in the Euclidean space with coordinates {x, t, z}, which we de- 
note as 7r X2 (Ci); (2) the curve defined in the xz-plane by gcd(lcoeff t (h), lcoeff f (/)) 
The equation of ir xz (Ci) is clearly h(t,z) = 0. Thus, by Proposition [T5l and 
Proposition [TBI and using elementary properties of the resultant, one gets that 
the considered values belonging to M.\ also satisfy h t (t, z) = 0, and hence they 
are contained in the set of real roots of Rest(/i, h t ); one may observe that this 
set contains also the real roots of gcd(lcoeff 4 (/i), lcoefff(/)). Arguing in a sim- 
ilar way for Ai 2 we would reach the condition Kes t (j,jt) = 0. Moreover, the 
(3)-values can be related with the asymptotes of the curves (in the tz-plane) 
h(t, z) = and j(t, z) = 0. So, the following theorem holds. 

Theorem 17 The (l)-values, (2)-values not corresponding to self-intersections 
of M., and (S)-values not corresponding to asymptotes of ii xz (C cr it) , are among 
the finitely many real roots ofRes t (h, h t ), Res t (j,j t ). 



5 Conclusions 

In this paper we have presented an algorithm for computing a critical set 
of a family of rational curves depending on a parameter. From the critical 
set, the topology types in the family can be derived. The algorithm is based 
on a geometric interpretation of known results for the implicit case, and on 
advantages of parametric representation, and requires certain properties on 
the family to be analyzed. These properties can be algorithmically checked. 
In our experimentation, we have found that the timings of the parametric 
algorithm are usually quite better than those of the implicit algorithm; in fact, 
the parametric algorithm is able to manage inputs that the implicit algorithm 
cannot deal with. On the other hand, the drawback of the provided algorithm 
is that it may determine critical sets bigger than those determined by the 
implicit algorithm, therefore containing superfluous values with respect to the 
implicit critical set. So, as a potential future line of research, one could address 
the problem of reducing the size of the output, trying to approach optimality. 
Furthermore, the method applies with exact coefficients. So, it would also 
be nice to consider the (challenging) possibility of applying it in the case of 
approximate coefficients. 

Acknowledgements. The author wants to thank the referees of the paper for 
their useful observations and suggestions, that helped to improve the paper. 
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6 Appendix I: Parametrizations of the families used in the com- 
parison table. 

Family 1: 

u := -78t 4 + 62t 3 + lit 2 + 88t + 1 + A(30t 4 + 81t 3 - 5t 2 - 28t + 4) 
v : = -lit 4 + 10t 3 + 57t 2 - 82t - 48 + A(— lit 4 + 38t 3 - 7t 2 + 58t - 94) 

The implicit equation has degree 4 (as a polynomial in x, y) and 95 terms; the 
total degree can be found (for this family and for the rest of the families in 
this appendix) in the comparison table provided in Subsection HI 

Family 2: 

u ■= 50 - 85t 5 - 55t 4 - 37t 3 - 35t 2 + 97t + A(-59 + 79t 5 + 56t 4 + 49t 3 + 63t 2 + 57t) 
v ■= -62 + 45t 5 - 8t 4 - 93t 3 + 92t 2 + 43t + A(-61 + 77t 5 + 66t 4 + 54t 3 - 5t 2 + 99t) 

The implicit equation has degree 5 (as a polynomial in x, y) and 161 terms. 
Family 3: 

2t 7 + t 5 A + A 2 t 4 - 3t 3 A 2 + 3A 3 t 2 - t 3 + 3tA - 2A 2 - 2t 4 A + t 6 A - t 4 
(t 4 -2tA + A 2 )(t 3 -tA + A 2 ) 

t 3 + At 2 - 1 
V : ~ t 3 - tX + A 2 

The implicit equation has degree 7 (as a polynomial in x, y) and 343 terms. 
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Family 4: 



_ -35 - 85* 3 - 55t 2 - 37t 66 + 43t 3 - 62t 2 + 77t 

U '~ 56 + 97t 3 + 50t 2 + 79t + -61 + 54t 3 - 5t 2 + 99t 
31 - 50t 3 - 12t 2 - 18* -59 + 49* 3 + 63* 2 + 57* 

V '~ 56 + 97* 3 + 50* 2 + 79* + -61 + 54* 3 - 5* 2 + 99* 

The implicit equation has degree 6 (as a polynomial in x, y) and 84 terms. 
Family 5: 

u ■= 50 - 85* 5 - 55* 4 - 37* 3 - 35* 2 + 97* + A(-59 + 79* 5 + 56* 4 + 49* 3 + 63* 2 + 57*) 
v ■= -62 + 45* 5 - 8* 4 - 93* 3 + 92* 2 + 43* + A(-61 + 77* 5 + 66* 4 + 54* 3 - 5* 2 + 99*) 



The implicit equation has degree 5 (as a polynomial in x, y) and 161 terms. 
Family 6: 



u ■= -47* - 91A 2 - 47* 3 - 61A 4 + 41* 5 - 58* 2 A 3 



v 



23* 2 - 84* 3 A + 19* 2 A 2 - 50*A 3 + 88* 5 A - 53* 2 A 4 



The implicit equation has degree 5 (as a polynomial in x, y) and 191 terms. 
Family 7: 

* 5 - t 2 X 2 - t - 2A + 1 

7/ " = 

t 3 -t 2 + t\- A 2 
* 5 + * 2 A-*-2A 2 + 1 

V : ~ * 3 - * 2 + *A - A 2 



The implicit equation has degree 5 (as a polynomial in x, y) and 219 terms. 
Family 8: 

-7t + 58* 2 - 94*A - 68* 3 + 14* 2 A - 35A 3 

-14 - 9t - 51A - 73* 2 - 73*A - 91A 2 
-50 + 50A + 67* 2 - 39*A + 8A 2 - 49*A 2 

V '~ -14 - 9* - 51A - 73* 2 - 73*A - 91A 2 



The implicit equation has degree 3 (as a polynomial in x, y) and 72 terms. 
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Family 9: 



-5 + 99t - 61A - 50A 3 - 12t 6 - 18A 6 
31 - 2Qt - 62A + t 2 - A7tX - 91A 2 
v := -1 + 94t 2 + 83A 2 - 86tA 2 + 23A 3 - 84t 3 A 

The implicit equation has degree 6 (as a polynomial in x, y) and 204 terms. 
Family 10: 

-85 - 55t - 37A - 35t 2 + 97t\ + 50A 2 

79 + 56t + 49A + 63t 2 + 57tA - 59A 2 
45 - St - 93A + 92t 2 + 43tA - 62A 2 

V '~ 79 + 56t + 49A + 63t 2 + 57tA - 59A 2 

The implicit equation has degree 2 (as a polynomial in x, y) and 30 terms. 
Family 11: 

97A + 50tA + 79A 2 + 56t 3 + 49tA 2 + 63A 3 

—93* + 92A + 43tA - 62t 3 + 77t A 2 + 66A 3 
-12 - 18t + 31A - 26tA - 62A 2 + i 2 A 

V '~ -93t + 92A + 43tA - 62t 3 + 77tA 2 + 66A 3 

The implicit equation has degree 3 (as a polynomial in x, y) and 97 terms. 
Family 12: 

57t - 59tA + 45A 2 - 8t 3 - 93t\ 2 + 92t 2 A 2 

' -18t + 3H 2 -26tA-62t 3 + t 2 A 2 -47A 4 
-1 + 94i 2 + 83A 2 - 86tA 2 + 23A 3 - 84t 3 A 

V '~ -18t + 31t 2 - 26iA - Q2t 3 + t 2 A 2 - 47A 4 

The implicit equation has degree 3 (as a polynomial in x, y) and 146 terms. 
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